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and which we shall refer to as multiple polylogarithms. When fc = 0, we define 
Q>^ I X{{}) := 1, where {} denotes the empty string. When k = 1, note that 
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Abstract. Historically, the polylogarithm has attracted speciaUsts and non- 
nJ ' specialists alike with its lovely evaluations. Much the same can be said for 

C ^ , Euler sums (or multiple harmonic sums), which, within the past decade, have 

arisen in combinatorics, knot theory and high-energy physics. More recently, 

we have been forced to consider multidimensional extensions encompassing the 

classical polylogarithm, Euler sums, and the Riemann zeta function. Here, 

^^ ^ we provide a general framework within which previously isolated results can 

^^ , now be properly understood. Applying the theory developed herein, we prove 

C J ' several previously conjectured evaluations, including an intriguing conjecture 

• ' of Don Zagier. 

fi: 

1. Introduction 
We are going to study a class of multiply nested sums of the form 

k ^^ ' / \ oo k / k 












is the usual polylogarithm E^ Q when s is a positive integer and IH ^ 1- Of 



course, the polylogarithm (L2) reduces to the Riemann zeta function 



(1.3) ^(,) = ^1, sR(s)>l, 

when 6=1. More generally, for any fc > the substitution Uj — X]i=7 ^i shows 



that our multiple polylogarithm (1.1) is related to Goncharov's [p5| by the equation 

Lisfc....,si(a:fe,--- ,a;i) = A i''"' j^ where yj:=TTa:^r\ 

\yi,---,ykj fj: 
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and 

(1.4) 



Lis,,...,si(a:^fe,.-. ,a;i) := ^ Y[ 

ni>->nfc>0 j=l 



"j ^J 



With each Xj = 1, these latter sums (sometimes called "Euler sums"), have been 
studied previously at various levels ol generality |, |6[ 0, g |l|, |l|, ^, |l|, |l|, |3|, |9[ 
E2I pU, M , the case k — 2 going back to Eule r pTf . Recently, Euler sums have arisen 
in combinatorics (analysis of quad-trees pu, [46[ and of lattice reduction algorithms 
|23[), knot theory ||l^, ^ |l^, |4^, and high-energy particle physics Q (quantum 
field theory). There is also quite sophisticated work relating polylogarithms and 
their generalizations to arithmetic and algebraic geometry, and to algebraic K- 
theory g 0, ^, ||, ||, ||, ^, ^, g. 

In view of these recent applications and the well-known fact that the classical 
polylogarithm (|l.2| ) often arises in physical problems via the multiple integrati on o f 
rational forms, one might expect that the more general multiple polylogarithm (1.1) 
would likewise find application in a wide variety of physical contexts. Nevertheless, 
lest it be suspected that the authors have embarked on a program of generaliza- 
tion for its own sake, let the reader be assur ed that it was only with the greatest 
reluctance that we arrived at the definition (LI). On the one hand, the polyloga- 
rithm ( |l.2| ) has traditionally been studied as a function of b with the positive integer 
s fixed; while on the other hand, the study of Euler sums has almost exclusively 
focused on specializations of the nested sum ( |l.4| ) in which each Xj — ±1. However, 
we have found, in the course of our investigations, that a great deal of insight is 
lost by ignoring the interplay between these related sums when both sequences of 
parameters are permitted to vary. Indeed, it is our view that it is impossible to 



fully understand the sums (1.2-1.4) without viewing them as members of a broader 
class of multiple polylogarithms. 



That said, one might legitimately ask why we chose to adopt the notation (1.1) 
in favour of Goncharov's (1.4), inasmuch as the latter is a direct generalization of 
the Li„ notation for the classical polylogarithm. As a matter of fact, the nota- 
tion ( |l.4[ ) (with argument list reversed) was our original choice. However, as we 
reluctantly discovered, it turns out that the notation (1.1), in which the second 
row of parameters comprises the reciprocated running product of the argument 



list in (1.4), is mo re s uitable for our purposes here. In particular, our "running 



prod uct" notation (1J-), in addition to simplifying the iterated integral representa- 
tion (|4.9| ) ( cf. pdj Theorem 16) and the various duality formulae (Section g — see eg. 
equ ation s (|6.7| ) and ( p.q )), brings out much more clearly the relationship (Subsec- 



tion g_^ between the partition integral (Subection 4.1), in which running products 
necessarily arise in the integrand; and "stuffles" (Subsections 5J, 5^). It seems also 
that boundary cases of certain formulae for alternating sums must be treated sepa- 
rately unless running product notation is used. Theorem p.5| with n = (Section g|) 
provides an example of this. 

Don Zagier (see eg. ^^) has argued persuasively in favour of studying special 
values of zeta functions at integer arguments, as these values "often seem to dic- 
tate the most important properties of the objects to which the zeta functions are 
associated." It seems appropriate, therefore, to focus on the values the multiple 
polylogarithms (^) take when the Sj are restricted to the set of positive integers, 
despite the fact that the sums (1.1) and their special cases have a rich structure as 
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analytic functions of the complex variables Sj. However, we allow the parameters 
bj to take on complex values, with each \bj\ > 1 and (61, si) ^ (1,1) to ensure 
convergence. 

Their importance notwithstanding, we feel obliged to confess that our interest in 
special values extends beyond mere utilitarian concerns. Lewin ^£| (p. xi) writes 
of a "school-boy fascination" with certain numerical results, an attitude which we 
whole-heartedly share. In the hope that the reader might also be convinced of the 
intrinsic beauty of the subject, we offer two modest examples. The first pq, M, 

00 k _. 2k 

vi,... ,Vk = l ] = 1 ^ ■' 

generalizes Euler's celebrated result 



°° 1 2 



.=1" 



and is extended to all even positive integer arguments in 0. The second (see 
Corollary ^ of Section ^), 



00 k , ^\ui + l °° ^ 

•••,i^fc=i j=i -^ vi,...,v^=i J 

(log 2)'= 



iy, + --- + Uk ,^^.^^^1 f^\ 2-^ (;.,+••■ + i^k) 



< fc e Z 



kl 

can be viewed as a multidimensional extension of the elementary "dual" Maclaurin 
series evaluations 






1/2" 

and leads to deeper questions of duality (Section o) and computational issues related 
to series acceleration (Section ^ . We state additional results in the next section and 
outline connections to combinatorics and g-series. In Section 0, we develop several 
different integral representations, which are then used in subsequent sections to 
classify various types of identities that multiple polylogarithms satisfy. Sections H 
through O conclude the paper with proofs of previously conjectured evaluations, 
including an intriguing conjecture of Zagier [ |69| and its generalization. 

2. Definitions and Additional Examples 
A useful specialization of the general multiple polylogarithm (^j), which is at 



the same time an extension of the polylogarithm (1^), is the case in which each 
bj = b. Under these circumstances, we write 



(2.1) M.„...,.,):=Ar;"--'^M= E n^-n^ 

and distinguish the cases b = 1 and b — 2 with special symbols: 
(2.2) C:=Ai and 6:^X2. 
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The latter (5- function represents an iterated sum extension of the poly logarithm (1.2) 
with argument one-half, and will play a crucial role in computational issues (Sec- 
tion]^) and "duality" identities such as (p. The former coincides with (1.4) when 



fc > 0, each Xj — 1, and the order of the argument list is reversed, and hence can 
be viewed as a multidimensional unsigned Euler sum. We will follow Zagier [ p9| 
in referring to these as "multiple zeta values" or "MZVs" for short. By specifying 



each bj = ±1 in (1.1), alternating Euler sums are recovered, and in this case, it 
is convenient to combine the strings of exponents and signs into a single string with 
Sj in the jth position when bj = +1, and Sj— in the jth position when bj = —1. 
To avoid confusion, it should be also noted that in H] the alternating Euler sums 
were studied using the notation 

ni>--->rifc>0 j = l 

where Si, . . . , s^ are non-zero integers and aj :— signum(sj). 

Additionally, n repetitions of a substring U will be denoted by [/". Thus, for 
example, 

A(P-,in:=A(f; \\)^ t n '-""""' 



(Et2j-1 ^^) (Et2j ^« 



Unit Euler sums, that is those sums (1.1) in which each Sj — 1, are also important 
enough to be given a distinctive notation. Accordingly, we define 

/I 1 \ °° '^ / '' 

(2.3) M^i.---^^^) -A '■••'= E nrME-' 

' }^l,...Mk=l J = l ^ l=j 

To entice the reader, we offer a small but representative sample of evaluations 
below. 

Example 2.1. Euler showed that 

oo _. n— 1 -. oc -. 

n—l fc— 1 n—1 

and more generally ||2^, |5^, that 

m-2 

2C(m, 1) = mC{m + 1) - E "^("^ ^ ^'^''^^ + -^)' 2 < m G Z. 

k=l 

The continued interest in Euler sums is evidenced by the fact that a recent American 
Mathematical Monthly problem [2^ effectively asks for the proof of C(2, 1) — C(3). 

Two examples of non-alternating, arbitrary depth evaluations for all nonnegative 
integers n are provided by 

Example 2.2. 

27r^" 

c({3,in-4-"c({4n = 



(4n + 2)!' 
previously conjectured by Don Zagier 169] and proved herein (see Section O); and 
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Example 2.3. 



n . 

C(2, {1, 3}") = 4-" 5](-l)'=C({4}"-^-) (4fc + l)C(4fc + 2) 

k 

-4^C(4j-l)C(4fc-4j + 3)j>, 
i=i 
conjectured in M and proved by Bowman and Bradley pi[. 



Example 2.4. An intriguing two-parameter, arbitrary depth evaluation involving 
alternations, conjectured in m and proved herein (see Section ||), is 



M({-1}"\ 1, {-!}") - (-1)™+^ E ^fc+n+l^™-fe 



(2.4) 



A;=0 



n + fc 



fc=0 ^ 



where 

(2.5) ^,:=Li.(i)=<5(r)=^-— , P, 



(log 2)'- 



k=l 



, Z, := (-l)'-C(O- 



The formula (2.4) is valid for all nonnegative integers m and n if the divergent 
m — case is interpreted appropriately. 



Example 2.5. If the Sj are all nonpositive integers, then 
Consequently, 

/ \ cxo A; /A; 

j 



Ui—0 



(2.6) 






j=i 1^ &j exp [J2i=i u^) 



- 1 



In particular, ( p.6[ ) implies 

(2.7) 



A 



0,...,0 
bi,...,b, 



/ .7 = 1 ^ 



Despite its utter simplicity, (2.7) points the way to deeper waters. For example, if 
we put bj — q^^ for each j = 1,2, ... ,k and note that 

^^ '^ ' '^ / ni>ri2>--->rifc>0 j = l 
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then ( |2.7D implies the generating function equaUty 

fe=o \y y ' 'y / „^i j,^g j^i y 

which experts in the field of basic hypergeometric series will recognize as a q- 
analogue of the exponential function and a special case of the g-binoniial theorem, 
usually expressed in the more familiar form 132] as 

°o fc(fc+l)/2 

{-zq; q)oc = V —, ^ — z''. 



The case fc = 1, 61 = 2, si = -n of (^ yields the numbers ||63| (A000629) 

(2.8) J(_„)^A2(-n) = ^^ = Li_„(i), 0<neZ, 

fe=i 
which enumerate Eq| the combinations of a simplex lock having n buttons, and 
which satisfy the recurrence 

d{~n) = l + J2[-)s{-j), l<neZ. 

Also, from the exponential generating function 

^ ^ " n! 2-e^ 2 - e^ 

n— 

we infer |g^, Q that for n > 1, ^S{—n) also counts 

• the number of ways of writing a sum on n indices; 

• the number of functions / : {1, 2, . . . , n} -^ {1, 2, . . . , n} such that if j is in 
the range of /, then so is each value less than or equal to j; 

• the number of asymmetric generalized weak orders on {1, 2, . . . , n}; 

• the number of ordered partitions (preferential arrangements) of {1, 2, . . . , n}. 

The numbers ^S{—n) also arise [|j| in connection with certain constants related 
to the Laurent coefficients of the Riemann zeta function. See [^ (A000670) for 
additional references. 

3. Reductions 



Given the multiple polylogarithm (IT), we define the depth to be fc, and the 



weight to be s := si + ■ ■ • + s^. We would like to know which sums can be expressed 
in terms of lower depth sums. When a sum can be so expressed, we say it reduces. 
Especially interesting are the sums which completely reduce, i.e. can be expressed 
in terms of depth-1 sums. We say such sums evaluate. The concept of weight is 
significant, as all our reductions preserve it. More specifically, we'll see that all our 
reductions take the form of a polynomial expression which is homogeneous with 
respect to weight. 

There are certain sums which evidently cannot be expressed (polynomially) in 
terms of lower depth sums. Such sums are called "irreducible" . Proving irre- 
ducibility is currently beyond the reach of number theory. For example, proving 
the irrationality of expressions like C(5, 3)/C(5)C(3) or C(5)/C(2)C(3) seems to be 
impossible with current techniques. 
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3.1. Examples of Reductions at Specific Depths. The functional equation 
(an example of a "stuffle" - see Sections ^ through |5^ ) 

reduces ({s, s). 

One of us (Broadhurst), using high- precision arithmetic and integer relations 
finding algorithms, has found many conjectured reductions. One example is 

(3.1) C(4, 1, 3) - -C(5, 3) + |iC(8) - |C(5)C(3) + K(3)'C(2), 

which expresses a multiple zeta value of depth three and weight eight in terms of 
lower depth MZVs, and which was subsequently proved. Observe that the combined 



weight of each term in the reduction (3.1) is preserved. The easiest proof of (|3.1 



uses Minh and Pctitot's basis of order eight | |55{ . 

Broadhurst also noted that although C(4, 2, 4, 2) is apparently irreducible in terms 
of lower depth MZVs, we have the conjectured|^ weight-12 reduction 

C(4, 2,4,2) = -ifiA(9-,3) - '-^C^ - ^({^,3) - f C(9)C(3) 
(3.2) - 1|0C(7)C(5) + 2C(6)C(3)2 + 14C(5, 3)C(4) 

+ 70C(5)C(4)C(3)-iC(3)' 

in terms of lower depth MZVs and the alternating Euler sum A(9— ,3). Thus, 
alternating Euler sums enter quite naturally into the analysis. And once the alter- 
nating sums are admitted, we shall see that more general polylogarithmic sums are 
required. 

We remark that the depth-two sums in (|3.2|), namely A(9— ,3), ^(9, 3), and 
^(5, 3), are almost certainly irreducible. For example, if there are integers ci, C2, C3, 
Ci (not all equal to 0) such that CiC(5, 3)-hC27rS-hC3C(3)^C(2) + C4C(5)C(3) = 0, then 
the Euclidean norm of the vector (01,02,03,04) is greater than 10^°. This result 
can be proved computationally in a mere 0.2 seconds on a DEC Alpha workstation 
using D. Bailey's fast implementation of the integer relation algorithm PSLQ p9[ , 
once we know the four input values at the precision of 200 decimal digits. Such 
evaluation poses no obstacle to our fast method of evaluating polylogs using the 
Holder convolution (see Section |^) . 

3.2. An Arbitrary Depth Reduction. In contrast to the specific numerical re- 
sults provided by ( |3.1[ ) and ( [3.2[ ) , reducibility results for arbitrary sets of arguments 
can be obtained if one is prepared to consider certain specific combinations of MZVs. 
The following result is typical in this respect. It states that, depending on the par- 
ity of the depth, either the sum or the difference of an MZV with its reversed-string 
counterpart always reduces. Additional reductions, such as those alluded to in 
Sections 1 and 2, must await the development of the theory provided in Sections 
4-7. 

Theorem 3.1. Let k be a positive integer and let si, S2, ■ ■ ■ , Sk be positive integers 
with si and s^ greater than 1. Then the expression 

C(S1, S2, • ■ • , Sk) + (-l)''C(Sfe, • ■ • , S2, Si) 

reduces to lower depth MZVs. 



^Both sides of (3.2) agree to at least 7900 significant figures. 
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Remark 3.2. The condition on si and Sk is imposed only to ensure convergence of 
the requisite sums. 

k 

Proof. Let N :— (Z+) denote the Cartesian product of k copies of the positive 
integers. Define an additive weight-function w : 2^ — » R by 



KA) := E n 



rJSA j=l 

where the sum is over ah n = {ni,n2, . . . , Uk) ^ A C N. For each I < j < k — 1, 
define the subset Pj of N by 

Pj :~ {fi (z N : rij < nj+i}. 

The Inclusion-Exclusion Principle states that 

(3.3) u;(nA^\p,)= Y. (-i)'^'^(n^. 

^i=i ^ TC{i,2,...,fc-i} ^jeT 

We remark that the term on the right-hand side of (^^) arising from the subset 
T = {} is C('5i)C(s2) ■ ■ ■ C(sfe) by the usual convention for intersection over an empty 
set. Next, note that the left-hand side of ( |3.3| ) is simply C(si, S2, . . . , Sk)- Finally, 
observe that all terms on the right-hand side of ( |3.3| ) have depth strictly less than 
k — except when T = {1, 2, . . . , A: — 1}, which gives 

k 

(_l)fc-i Y^ Y^ nj'' = (-l)'=-iC(sfc, . . . , S2, si) + lower depth MZVs. 

"i<"2<---<"fc i=i 



This latter observation completes the proof of Theorem 3.1. D 



4. Integral Representations 
Writing the definition of the gamma function [B9| in the form 

/•oc 

r-'T(s)=/ {logxy-'^x-'-^dx, r > 0, s > 0, 



Ji 
it follows that if each Sj > and each \bj\ > 1, then 

1 /-"^ (logx)^!-!^/ S2,...,Sfc \dx 



r(si) Ji a;6i - 1 \^^2, • • • , a;6fc/ a; ' 

a representation vaguely remindful of the integral recurrence for the polylogarithm. 
Repeated application of (4.1) yields the /c-dimensional integral representation 



<") M ;•■■' = ■■■ n 



k 



si,...,sA r r-TT (logxjp-^dx 



6i, . . . , bkj Ji J I f}^ n^Mbj riLi ^r - i)^j ' 
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which generahzes Crandall's integral EQ] for C(si, • ■ ■ , Sfc). An equivalent formula- 
tion of (U) is 

(s^ si\ f^ r°° -X u''~^du, 



bi,... ,bkj Jo Jo fj, T{s, ) (6j exp ( ELi "0 - 1) 



the integral transforms in ( |4.3| ) replacing the der ivat ives i n (|2.6| ) . 

Although depth- dimensional integrals such as (|4.2| ) and ( [4. 3D are attractive, they 
are not particularly useful. As mentioned previously, we are interested in reducing 
the depth whenever this is possible. However, since the weight is an invariant of all 
known reductions, we seek integral representations which respect weight invariance. 
As we next show, this can be accomplished by selectively removing logarithms from 
the integrand of (12), at the exp ense of increasing the number of integrations. At 
the extreme, the representation (4.2) is replaced by a weight- dimensional integral 
of a rational function. 

4.1. The Partition Integral. We begin with the parameters in ([L1|). Let Ri, 
i?2, . . . , Rn be a (disjoint) set partition of {1,2,..., k}. Put 

If di, ^2, . . . , d„ are real numbers satisfying \dm\ > 1 for all m and ridi ^ 1, then 

ri,...,r 
di , . . . , d„ 



Vl,... ,l^n — l 7n—l ^ j — 7n 



t^l ,... ,t/„ — 1 771—1 i^Rrix j—m 

i: no-/ -J n .!'°si....t . 

1^1,... ,i/„ = l m=l "^1 "^1 iefl™ 



r(s,)xi~ 



Now collect bases with like exponents and note that "nm=i Iliei? — 117=1 •" ^^ 
follows that 



m:::::;:::)^/ i \ ^ n<-.-nn 



^ A (loga:^^)'''' ^dx 

X 



(4.4) 



1^1,... ,1/^ — 1 m— 1 j — li^Hj 

}][ r(s,)x, 

1 '^ m=l ^ j = l ieflj 
^ JT (^OgX^Y'^^ dXj 



j.^1 r(sj)xj 

on summing the n geometric series. 

Example 4.1. Taking n — k, we have Rm = {"ti}, and r„j = s,„ for all 1 < ?7i < n. 
In this case, (4.4) reduces to the depth-dimensional integral representation (4.2). 



10 
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Example 4.2. Taking n = 1, we have Ri = {1,2,..., k} and ri 



EL^.-if 



we also put d :— Yii^i ^ji then (14) yields the seemingly wasteful fc-dimensional 
integral 



'0^'&. 



for a polylogarithni of depth one. 



. .n^-O^n^^^^iii^ 



Example 4.3. Let Sj — 1 for each 1 < j < fc, ro = and let ri,r2, . . . ,r„ be 
arbitrary positive integers with X]m=i ^™ ^ ^- ^'^^ 1 < m < n define 






In this case, (iA) yields a weight-dimensional integral of a rational function in k 
variables: 



(4.5) 



A 



n 



1 ■ ■ ■ J ' n 



ti, . . . , fij^ 



nH™n--i n 



4=1 



da;. 



J = l ■' 



where u„i = X]i=i ^j- '^'^ interesting specialization of (4.5) is 

dx dy dz 

^^^='i k 



C(2,l) = 



CXD /"OO /"OO 



C(3). 



1 xyz{xy-l){xyz-l) 



OO /"OO /'OO 



dx dy dz 
xyz{xyz — 1) 



Although it may seem wasteful, as in Example |4.1| above, to use more integra- 
tions than are required, nevertheless such a technique allows an easy comparison 
of multiple polylogarithms having a common weight but possessing widely differing 
depths. For example, from the four equations 

^fs + t\_ 1 [^ f°" (iogxy-^ilogyf-^dxdy 



(4.6) 



M'Ul 



1 

1 



1 {abxy - l)xy 

°° '"^ (\ogxy-\\ogyY-^dxdy 



ab 
s, t 
a, abj r(s)r(t) Ji 

t,s\_ 1 _ 

b.abj r(s)T(t).l, A {by - l){abxy - l)xy ' 

°° '"^ {\ogxy~H^ogyY'^dxdy 
1 



I {ax — l){abxy ~ l)xy 
°° '"^ {\ogxy-H^ogyY~^dxdy 



r(s)r(t) J, 
1 



yaj \bj r(s)r(t) Ji 

and the rational function identity 

1 1 



{ax- l)(6y- l)xy 



(4.7) 



{ax — l){by — 1) abxy — l\ax — 1 by — 1 



1 



1 



1 



the "stufRe" identity (see Section 5.1) 

(4.8) 

follows immediately. The connection between "stufHe" identities and rational func- 



,a/ \bl \a,abl \b,abl \ ab 



tions will be explained and explored more fully in Section 5.3 
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4.2. The Iterated Integral. A second approach to removing the logarithms from 
the depth-dimensional integral representation ( |4.2| ) yields a weight-dimensional it- 
erated integral. The advantage here is that the rational function comprising the 
integrand is particularly simple. 

We use the notation of Kassel |Q for iterated integrals. For j — 1,2, . . . ,n, let 
fj : [a,c] ^ R and flj := fjijjj) dyj. Then 



J a A — ^^a 



Iahiyi)C^^---^ndy, ifn>0 
1 if n = 0. 



For each real number b, define a differential l-form 

dx 



With this definition, the change of variable y i— > 1 — y generates an involution 
w(6) I— > w(l — fe). By repeated application of the self-evident representation 

b"'m-'' = / w^-iy™-! dy, 1 < m G Z 

Jo 

one derives from (l]^) that 



i/i,... .Uk=ij=i 

k 

n 






(4.9) = {^If ff[ur'u{b,). 

Letting s := si + S2 + • • • + Sfc denote the weight, one observes that the represen- 
tation ( |4.9| ) is an s-dimensional iterated integral over the simplex 1 > yi > y2 > 
• • • > j/s > 0. Scaling by q at each level yields the following version of the linear 
change of variable formula for iterated integrals: 

("«) ^::;::::::) -<;:;:::;*> <-"'rn-»'-""<'^' 

for any real number g ^ 0. 

Having seen that every multiple polylogarithm can be represented ( |4.9| ) by a 
weight-dimensional iterated integral, it is natural to ask whether the converse holds. 
In fact, any convergent iterated integral of the form 

(4.11) / n^^M 

•^0 r=l 
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can always (by collecting adjacent wq factors - note that for convergence, a{s) ^ 0) 
be written in the form 

r-l fe , 



(4.12) / f[co:r'^ib,) = i-irx 



j^i M^---^bk 



where 



3 

(4.13) Q^bj=a[ ^s, 

j=i 



We remark that the iterated integral representation ([l.9|) and the weight-dimen- 



sional non-iterated integral representation (45) of Example O are equivalent under 



the change of variable Xj =yj-i/yj, 2/o := 1, j = 1, 2, . . . , s. In fact, every integral 



representation of Section [4.1| has a corresponding iterated integral representation 



under the aforementioned transformation. For example, the depth-dimensional 



integral (4^) becomes 



^/si,...,Sfe\ yr fy^ ^ {log{yj-i/yj)Y' ^ dyj 



bi,...,bkj fj^ Jo r(sj ) {bj ~yj) 

The explicit observation that MZVs are values of iterated integrals is apparently 
due to Maxim Kontscvich [^. Less formally, such representations go as far back 
as Euler. 



5. Shuffles and Stuffles 

Although it is natural to study multiple polylogarithmic sums as analytic objects, 
a good deal can be learned from the combinatorics of how they behave with respect 
to their argument strings. 

5.1. The Stuffle Algebra. Given two argument strings s = {si, . . . ,Sk) and t = 
{ti, . . . ,tr), we define the set stufile(s, i) as the smallest set of strings over the 
alphabet 

{Si,.. .,Sk,ti,.. . ,tr, + , , , [ , ) } 

satisfying 

• (si,. . . ,Sfc,ti, . . . ,tr) e stuffle(s, t). 

• If a string of the form ([/, s„, i™, V) is in stuffle(s, i), then so are the strings 
{U,tm,Sn,V) and (t/, S„ + tm,V). 

Let a ~ (ai, . . . , Uk) and b = {bi, . . . ,br) be two strings of the same length as s and 
t, respectively. We now define 



(5.1) 5r:=5r,^, 

\a, b 

to be the set of all pairs (Jt) with u G stuffle(s, i) and c — (ci, C2, . . . , Ch) defined as 
follows: 

• h is the number of components of u, 

• Co := ao := bo := 1, 
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• for 1 < j < h, if Cj-i — an-ibm-i, then 

O^n^m—l: II ^j — Sn^ 
a-n-lbm, iiuj = tm- 

5.2. Stuffle Identities. A class of identities which we call "depth-length shujfles" 
or "stujfle identities" is generated by a formula for the product of two A-functions. 
Consider 

<:>G^^{ t n«7"'(i:".)""}{ £ n^-'(i:«. 

If we put 



n., : 



then it follows that 

mi > ■ ■ ■ > mj. > 

Rewriting the previous expression in terms of A-functions yields the stuffle formula 

where the sum is over all pairs of strings (^j G ST(t'i^. 
Example 5.1. 

\a, bj \cj \a, b, be) \ a, be 

\a, ac, be) \ ac, be J \c, ac, be 

When specialized to MZVs, this example produces the identity 

C(r, s)C(i) = C(^ s, t) + C(r, s + t)+ C(r, i, s) + C(r + t, s) + C{t, r, s). 

The term "stufHe" derives from the manner in which the two (upper) strings are 
combined. The relative order of the two strings is preserved (shuffles), but elements 
of the two strings may also be shoved together into a common slot (stuffing), thereby 
reducing the depth. 



5.3. Stuffles and Partition Integrals. In Section 4.1, an example was given 



in which a stuffle identity (fl^) was seen to arise from a corresponding rational 



function identity (4.7) and certain partition integral representations (|4.6|). This is 



by no means an isolated phenomenon. In fact, we shall show that every stuffle 



identity is a consequence of the partition integral (4^) applied to a corresponding 
rational function identity. 

Theorem 5.2. Every stujfle identity is equivalent to a rational function identity, 
via the partition integral. 
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Before proving Theorem 5.2, we define a class of rational functions, and prove 
they satisfy a certain rational function identity. Let s = (si,...,Sfc) and t — 
(ii,...,tr) be vectors of positive integers, a nd l et a — (ai,...,afe) and /? = 
(/3i, . . . ,Pr) be vectors of real numbers. As in (|5.l|), put 



ST^ST 



s.t 



d, (3 J 



and define 



T = T(a,/5):=^7: [ ^^] ^ ST 



Let / : r ^ Q [71 , 72 , ■ • ■ ] be defined by 



(5.3) 



f{ii,---nh):=X{{i,-i)-'. 



Then we have the following lemma. 



Lemma 5.3. Let f be defined as in (5.S). Then 



/(«)/(/3)= E /(^)- 



Proof of Lemma 5.c. Apply ( p.2| ) with a — a and b — (3. In view of (2.7), the 
lemma follows on taking s and t to be zero vectors of the appropriate lengths. 



Proof of Theorem 5.i. Let s, i, a, and b be as in (5.2). Let a and (5 be given by 



aj:=aj\Yxu Pj.= bj'Wy,. 
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Applying Lemma 5.3 and the partition integral representation (4.4) to the depth- 
dimensional integral (4.2) yields 



m:m H -L /wn 



M _ / P r ff-\ rr C^ogx jY^-^dx 



a/ \bj [Ji Ji jj^ ^(.Sj)xj 



oo ^ n^„n,M,-i 



fwU 



(logyj)*^ ^dyj 



,.^1 nwvj 



r 

n 



(logyjY^-^dyj 



j^i ^(.tj)yj 



.s.r-r«Mn^^^fS?^} 



7eT(5,/3) i=i 

-TT (l0gyj)*^~^^yj 

/i ntj)yj 



(?)e^r(;:s) 



as required. 



5.4. The Shuffle Algebr a. A s opposed to depth-length shuffles, or stuffles, which 
arise from the definition (lA) in terms of sums, the iterated integral representa- 
tion (4.£) gives rise to what are called em "weight-length shuffles", or simply "shuf- 
fles" . Weight-length shuffles take the form 



(5.4) 



Jo Jo Jo 



where the sum is over all (" „") permutations a of the set {1, 2, . . . , n + m} which 
satisfy a~^{i) < cr~^(j) for all 1 < « < j < n and n + 1 < i < j < n + m. In 
other words, the sum is over all {n + 77i)-dimensional iterated integrals in which the 
relative orders of the two strings of 1-forms ili, . . . , 57„ and r2„+i, . . . , fln+m are 
preserved. 



Example 5.4. 




C(2,1)C(2) 


Jq Jo 




= — 6 / WgWj* ~ 3 / ujQLJiUJoUjf — / ujQUjfuJouJi 
Jo Jq Jq 




= 6C(3,l,l)-t-3C(2,2,l) + C(2,l,2). 



In contrast, the stuffle formula gives 

C(2, 1)C(2) = 2C(2, 2, 1) + C(4, 1) -I- C(2, 3) -I- C(2, 1, 2). 

Note that weight-length shuffles preserve both depth and weight. In other words, 
the depth (weight) of each term which occurs in the sum over shuffles is equal to 
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the combined depth (weight) of the two muhiple polylogarithms comprising the 
product. 

Though it may appear that the shuffles form a rather trivial class of identities sat- 
isfied by iterated integrals, it is worth mentioning that the second proof of Zagier's 



conjecture (see Corollary H of Section 11.2) uses little more than the combinatorial 



properties of shuffles H . In addition, both shuffles and stuffles have featured in the 
investigations of other authors in related contexts |^, ^, |4^, ^, ^, |4[ |5^, ^, ^ 
MM. 



6. Duality 

In ^^, Hoffman defines an involution on strings si, . . . , Sfe. The involution co- 
incides with a notion we refer to as duality. The duality principle states that 
two MZVs coincide whenever their argument strings are dual to each other, and 
(as noted by Zagier |£^) follows readily from the iterated integral representation. 
In |1^], Broadhurst generalized the notion of duality to include relations between 
iterated integrals involving the sixth root of unity; here we allow arbitrary complex 
values of bj. Thus, we find that the duality principle easily extends to multiple 
polylogarithms, and in this more general setting, has far-reaching implications. 

6.1. Duality for Multidimensional Polylogarithms. We begin with the iter- 



ated integral representation (49) of Section |4. 4 Reversing the order of the omegas 
and replacing each integration variable y by its complement 1 — y yields the dual 
iterated integral representation 

J — K 

where again s = si + ■ ■ ■ + Sk is the weight. 
Example 6.1. Using (0), (^) , and del]) , we have 



A (^^'_\) - £ c.(0) a;(l) u;{-l) = - ^' c^(2) c.(0) a;(l) = -A Q| ^j , 
which is to say that 

„2 Z^ fc ~ Z^ n2" ^ A:2 ' 

n— 1 k—1 n—1 k—1 

a result that would doubtless be difficult to prove by naive series manipulations 
alone. 



When bi = b2 = ■ ■ ■ = bk = b, the two dual iterated integral representations (4.9) 
iimplify as follows: 

s,,...,Sk) = {-If C n ^r'^{b) = (-1)^+^ /' n Lo{l - b)u;'r\ 
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A somewhat more symmetric version of (6.2) is 

1 m 

(-irAb(,si + 2,{ir,...,s„, + 2,{ir'") = {-ly n^o'^' "'^' 



^fc 



(6.3) = (-1)^/ n^i-tv 



1 1 

rj + l Sj + 1 




]=m 



where r :— ^ ■ r^ and, as usual, s :— J2j -^j- 

6.2. Duality for Unsigned Euler Sums. Taking 6 = 1 in (6.3), we deduce the 
MZV duality formula (cf. |4) p. 483) 

(6.4) C(si + 2, {ir\ . . . , s™ + 2, {I}'-") = C(r™ + 2, {1}^", . . . , ri + 2, {1}^0 

for multidimensional unsigned Euler sums, i.e. multiple zeta values (MZVs) . 

Example 6.2. MZV duality ^^ gives Euler's evaluation C(2, 1) = ((3), as well 
as the generalizations C({2, 1}") = C({3}"), and C(2, {1}") = C(" + 2), valid for all 
nonnegative integers n. 



In [|60[ a beautiful extension of MZV duality (6.4) is given, which also subsumes 
the so-called sum formula 

^ C{ni,n2,... ,nk) = C{N), 

nj>Sj I 

conjectured independently by C. Moen [Q and M. Schmidt [Q, and subsequently 
proved by A. Granville [^. We refer the reader to Dr. Ohno's article for details. 
The duality principle has an enticing converse, namely that two MZVs with 
distinct argument strings are equal only if the argument strings are dual to each 
other. Unfortunately, although the numerical (and symbolic) evidence in support 
of this converse statement is overwhelming, it still remains to be proved. In the 
case of self-dual strings, the conjectured converse of the duality principle implies 
that such a MZV can equal no other MZV; moreover we find that certain of these 
completely reduce, i.e. evaluate entirely in terms of (depth-one) Riemann zeta 
functions. 

Example 6.3. The following self-dual evaluation, previously conjectured by Don 
Zagier ||69|] 

27r^" 
C({3,l}") = 4-"C({4r)^ ^^^^^^, , 0<nGZ, 

is proved herein (see Section O) . 
Example 6.4. The evaluation 

n , 

C(2, {1, 3}") = 4-" ^(-1)'=C({4}"-') (4fc + 1)C(4A: + 2) 
fc=0 '^ 

k . 

-4^C(4j-l)C(4fc-4j + 3)L 0<neZ 

conjectured in m and recently proved by Bowman and Bradley |1^ is also self-dual. 
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Example 6.5. The self-dual two-parameter generalization of Example 6.3 

2(m -M)7r^('"+i)"+2m 



C({2r,{3,{2}™,l,{2}™r) = 



< 771, n G Z, 



(2(7n + l)(277+l))! ' 
remains to be proved. 

We conclude this section with the following result, since the special case p = 1 
has some bearing on the MZV duality formula (6.4). 

Theorem 6.6. Let \p\ > 1. The double generating function equality 



oc oc 
m— n— 



^7n-\-l 7i-\-l 



2/"+Up(777 + 2,{in-2Fl 



y,-x 



holds. 



Proof. By definition (2.1) of A 



•pj 



EE^ 

m=0 71—0 



^7n-\~l n-{-l 



7/"+Up(777-H2,{l}") 



1 



fc-1 



m=0 fe=l ^ j = l 



n 1 



m=0 fc=l 

°° (y)fe / X 



{y)k 



f^m+lf.\pk 



E Ky)k 
kip'' 
fc=l ^ 



(y)fc(-a:)fc 
^ klp''{l ~x)k 



E 



2F1 



J/, -a; 
l-x 



as claimed. 



D 



Remarks 6.7. In pj it was noted that the p = 1 case of Theorem 3.1 is equivalent 
to the m — \ case of MZV duality (6.4) via the invariance of 

T{l-x)V{l-y) 



2F1 



(6.5) 



y-,-x 

l-x 



1 = 



T{l~x~v) 



= exp|f;(x^-+y'=-(x + y)'=)^| 



with respect to the interchange of x and y. However, it appears that this observation 
can be traced back to Drinfeld ||2^ . In connection with his work on series of Lie 
brackets, Drinfeld encountered a scaled version of the exponential series above, 
and showed that the coefficients of the double generating function satisfy Cmn = 
Cnm and Cmo = com cvaluatcs to C,{ra + 2), up to a so-called Oppenheimer factor 
which we omit ([Q, p. 468). In our notation, this is essentially the statement that 

C(777 + 2,{l}")=C(n + 2,{ir). 

Note that Theorem ^.l| in conjunction with (6.5) shows that Q{ra + 2, {1}") com- 
pletely reduces (i.e. is expressible solely in terms of depth-1 Riemann zeta values) 
for all nonnegative integers m and n. In particular, the coefficient of x"^~^y'^ gives 
Euler's formula (Example 2.1); and taking the coefficient of x"^^^y^ provides a 
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much simpler derivation of Markett's formula | |51| for C("^j li l)j m > 2. Thus, the 
complete reducibility of <^{m + 2, {l}") is a simple consequence of the instance ( |6.5[) 
of Gauss's 2F1 hypergeometric summation theorem ||, ^ ^. Wenchang Chu [|19| 
has elaborated on this idea, applying additional hypergeometric summation the- 
orems to evaluate a wide variety of depth-2 sums, including nonlinear (cf. [ pl[ ) 
sums. 

It would be interesting to know if there is a generating function formulation of 
MZV duality at full strength (3.4). Presumably, it would involve an analogue of 
Drinfeld's associator in 2m non-commuting variables. 



Sfe) : 


-AX 


■,Sk\ 


00 k / ^ 

- E n^- E- 

vi,... ,Vk = l j = l ^ i=j 



6.3. Duality for Unit Euler Sums. Recall the (5-function was defined (2^) as 
the nested sum extension of the polylogarithm at one-half: 

(6.6) (S(si, 

Due to its geometric rate of convergence, (5-values can be computed to high preci- 
sion relatively quickly. On the other hand, the unit Euler /i-sums ( |2.3| ) converge 
extremely slowly when the bj all lie on the unit circle. In particular, the slow 
convergence of the unit (±1) argument /i-sums initially confounded our efforts to 
create a data-base of numerical evaluations from which to form viable conjectures. 
Nevertheless, there is a close relationship between the (5-sums and the /i-sums, as 
we shall presently see. 

Taking 6 = 2 in (|6.3D, we deduce the "delta-to-unit-mu" duality formula 



(6.7) <5(si + 2,{ir,...,s,„ + 2,{ir-) 

= (-i)'-+x{-ir'"+i, {iY-+\ . . . , {-ir+Mir^+i). 

Thus, every convergent unit (±1) argument /^-sum can be expressed as a (rapidly 
convergent) (5-sum. The converse follows from the more general, but less symmetric 
formula, arising from ( |6.2| ): 

(6.8) <5(si, . . . , sfe) = (-l)V(-l, {lY"-', ■•■,-!, {1}^^"')- 



Example 6.8. 



^(1) = E ^ = - Ml) = E ^— = -M-i), 

v—l u—1 

and more generally, for all nonnegative integers n, we have 

oc ^ 

(6.9) 5{n + 1) = ^ ^^^^ = Li„+i(i) = -/.(-I, {!}"). 

Example 6.9. For all nonnegative integers n, 

(6.10) 5{{ir) = (-l)V({-ir) = (log2)"/n!, 

(6.11) ^(2, {in = (-i)"+i/i({-ir+\i), 

and more generally, 

(^dl}™, 2, {1}") = (-l)"+"+V({-l}"+\ 1, {-1}™), < m, n e Z. 
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Example 6.10. 

5{l,n+l)=n{-l,{iy\-l), 0<nGZ, 



and in particular, remembering (2^,2^ 6^) that 6{r) — Lir(^), we have 



,5(1,0) = l-log2=l-J(l), 

,5(1,2) = IS{2)d{l)~'jS{3) + 4^S^l). 

Integer relation searches (see [|o) or [0 for details) have failed to find a similar 
formula for (5(1,4). However, 

2n-l 

2^(1, 2n - 1) = Y^ {-iy+^5{j)S{2n - j), 1 < n G Z. 

Also, 

Sil,-n) = Y[^ , l<neZ, 



where the S{—iy) are the simplex lock numbers ( |2.8| ) and the i?i/ are the Bernoulli 
numbers R]. More generally, if ni is a positive integer and 712,713, ... ,71^ are all 
nonnegative integers, then 

C r Tj 

5{s,-nr,. 
where 



-712, -ril) = I n £ '^(''j) f^(* - i^r - 1), S e C, 



Tj := 7ij + i/j^i + 1, A{vj) := -—-J rj] Br^_^^ , Vq := -1. 

7. The Holder Convolution 

Richard Crandall [BTI (see also |22J) describes a practical method for fast eval- 
uation of MZVs. Here, we develop an entirely different approach which is based 
on the fact that any multiple polylogarithm can be expressed as a convolution of 
rapidly convergent multiple polylogarithms. We have used such representations to 
compute otherwise slowly convergent alternating Euler sums and (unsigned) MZVs 
to precisions in the thousands of digits. Lest this strike the reader as perhaps an 
excessive exercise in recreational computation, consider that many of our results 
were discovered via exhaustive numerical searches [[7| for which even hundreds of 
digits of precision were insufRcient, depending on the type of relation sought [ |lO| . 

A publicly available implementation of our technique is briefly described in Sec- 
tion 7.2. There are also interesting theoretical considerations which we have only 



begun to explore. See equations (0)-(7^) below for a taste of what is possible 



7.1. Derivation and Examples. We have seen how multiple polylogarithms with 
unit arguments can be expressed in terms of rapidly convergent (5-sums. What if the 



arguments are not necessarily units? In the iterated integral representation (4.9) 
the domain 1 > yj > yj+i > in s = "^jSj variables splits into s -I- 1 parts. 
Each part is a product of regions 1 > yj > yj+i > 1/p for the first r variables, 
and 1/p > yj > yj+i > for the remaining s — r variables. Next, yj 1-^ 1 — yj 
replaces an integral of the former type by one of the latter type, with 1/p replaced 
by 1/q := 1 - 1/p. 
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Motivated by these observations, we consider the string of differential f-forms 
which occurs in the integrand of the iterated integral representation (|4.9|) and define 



ar 



0, otherwise. 



Then 



Si,. 



,Sfc 

,bk 



1 s 



(7.1) 



r=0 ^""^ j=r 



Uj{l — Uj) 



1/p 



n ^"^"j) 

j=r+l 



Thus, by means of ( 4.11 ), ( 4.12 ), and ( 4.13 ), we have expressed the general multiple 
polylogarithm as a convolution of Ap with Xq for any p, q suc h th at the Holder 
condition 1/p + 1/q = 1 is satisfied. For this reason, we refer to ( [7.l|) a s the Holder 
convolution. Note that the Holder convolution generalizes duality (|6J ) for multiple 
poly logarithms, as can be seen by letting p tend to infinity so that ( 4.10 ) Xp -^ 0, 
and q ^ 1. 

MZV Example. For any p > 0, q > with 1/p + 1/q = 1, 

C(2, 1,2, 1,1,1) = Ap(2,l,2,l,l,l) + Ap(l,l,2,l,l,l)A,(l) 

+Ap(l, 2, 1, 1, 1)A,(2) + Ap(2, 1, 1, 1)A,(3) 

+Ap(l, 1, 1, 1)A,(1, 3) + Ap(l, 1, 1)A,(2, 3) + Ap(l, 1)A,(3, 3) 

+Ap(l)A,(4,3) + A,(5,3) 

= C(5,3). 

The pattern should be clear. For 1 < j < m, define the concatenation products 



Cat{s, + 2,{ir-} = {sj+2,{ip,. 



1 



Cat{r, + 2,{l}^'}-{r,+2,{l}^ 



i=j 



,n + 2,{ir^}, 



and Qm+i '■= bo := {}. Then the Holder convolution for the general MZV case is 
given by 

C(a™) = E EAp(s,+2~i,{ip,a,+i)A,({ir,6,^i) 
j=l [ t=o 



(7.2) 



E Ap({l}^ a,+i)A,(r, + 2 - t., {ir^ 6,_i) + A,(5,„) 



abm). 



Of course, a,„ and b„i are the dual strings in the MZV duality formula (6.4). Since 
the sum s Xp converge geometrically, whereas MZV sums converge only polyno- 
mially, ( |7.2D provides an excellent method of computing general MZVs to high 
precision with the optimal parameter choice p = q = 2. For rapid computation of 
general multiple polylogarithms, it is simplest to use the Holder convolution ( [7.l[ ) 
directly, translating the iterated integrals into geometrically convergent sums on 
a case by case basis, using ( |4.9D . 
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Alternating Example. 



A(2,l-) = / ojiO)Lu{l)uji-l) 

JQ 



1/p j-l/q j-l/p 

uj{0)lu{1)uj{-1) - uj{l) w(l)w(-l) 

"'0 "'0 

1/q pl/p pl/q 

uj{0)lu{1) Lo{-l)- I LO{2)UJ{0)UJ{1) 

Jo Jo 



Ap(2, 1-) + Ap(l, l-)Ag(l) + Ap(l-)A,(2) - A, 



1,2 
2,1 



1,2 
2,1 



Although we could now work out the explicit form of the analogue to ([7^) in 
the alternating case, the resulting formula is too complicated in relation to its 
importance to justify including here. 

In addition to the impressive computational implications already outlined, the 
Holder convolution (7.1) gives new relationships between multiple polylogarithms, 
providing a path to understanding certain previously mysterious evaluations. For 
example, taking p = q — 2 shows that every MZV of weight s can be written 
as a weight-homogeneous convolution sum involving 2s (5-functions. Furthermore, 
employing the weight-length shuffle formula (|]j) to each product shows that every 
MZV of weight s is a sum of 2* (not necessarily distinct) (5-values, each of weight 
s, and each appearing with unit (-1-1) coefficient. In particular, this shows that the 
vector space of rational linear combinations of MZVs is spanned by the set of all 
(5-values. Thus, 



C(3) 



1/2 .1/2 .1/2 



1/2 



1/2 .1/2 

UJlLUl / iUl 

Jo 



UJqUJiUJi 



5(3) 



1/2 



{ll>i ■ ajQUJl + LOq ■ UJl ■ LUl + UJqLl>i ■ UJl) 





1/2 



{lOiLUi ■ LOl + LOl ■ LUl ■ LOl + LOl ■ ^I'^-'l) + (5(2, 1) 

= (5(3) + (5(1, 2) + (5(2, 1) + (5(2, 1) + (5(1, 1, 1) + (5(1, 1, 1) + (5(1, 1, 1) + (5(2, 1). 

Polylog Example. Applying (7.1) to C(n-f 2), withp = q = 2 provides a lovely 
closed form for (5(2, {1}"). Indeed, 



n+2 



(7.3) 



C(n + 2) = (5(2, {1}") + J2 Sir)Si{ir+'-n. 



r=\ 



The desired closed form follows after rearranging the previous equation (7.3) and 
applying the definition (a.E) and the result ( 6.10| ) in the form 6{r) = Lir(^) and 
'^({1}'^) = {^og2y/r\, respectively. 
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Example 7.1. Putting n = 1 in ( |7.3| ) gives [g 

Ti=l n=l j = l ■' 



In fact, formula ( |7.3| ) is non-trivial even when n = 0. Putting n = in ( |7.3| 
gives the classical evaluation of the dilogarithm at one-half: 

(7.5) 2Li2(i)=C(2)-(log2)2 i.e. Y. = ^^^^^ - \[\o^2f . 



=1 



2"n2 



Differentiation of (7.1) with respect to the parameter p provides another avenue 
of pursuit which has not yet been fully explored. We have used this approach to 
derive (5(0, {1}") ~ (5({1}"), but in fact, removing the initial zero is trivial from 
first principles. 

7.2. EZ Face. A fast program for evaluating MZVs (as well as arithmetic expres- 



sions containing them) based on the Holder convolution formula (7.2) has been 
developed at the CECMn, and is available for public use via the World Wide Web 
interface called "EZ Face" (an abbreviation for Euler Zetas intcrFace) at the URL 



http: //www. cecm. sf u. ca/projects/EZFace/ 



This publicly accessible interface currently allows one to evaluate the sums 



z(si,... ,Sfc) := ^ ]^n^- "' a 

ni>--->nfc j = l 

for non-zero integers si, . . . , s/c and Oj :— signum(sj), and 
zp(p,si,... ,Sk) := Yl ^""'11 



k 
ni>--->nfc j=l 



for real p > I and positive integers si, . . . , Sfc. The code for evaluating these sums 
was written in C, using routines from GMP, the GNU Multiprecision LibraryP^. Our 
implementation permits the precision of the evaluation to be set anywhere between 
10 and 100 digits. Progress is currently underway to extend the scope of sums that 
can be evaluated. The exact status of the EZ Face is at any moment documented 
at its "Definitions" and "Using EZ-Face" pages. 

In addition to the functions z and zp, the lindep function, based on the LLL 
algorithm |Q for discovering integer relations ||l^ satisfied by a vector of real num- 
bers, can be called. An integer relation for a vector of real numbers {xi, . . . , Xn) is a 
non-zero integer vector (ci, . . . , c„) such that J27=i '^i^i ~ 0- The required syntax is 
lindep([a;i, . . . , a;„]), where xi, ... , a;„ is the vector of values for which the relation 
is sought. One must ensure that the vector of real numbers is evaluated to sufficient 
precision to avoid bogus relations and other numerical artifacts. The lindep code 
was written by Michael Monagan and Greg Fee, both of the CECM, and is available 
on request. Send e-mail to either monagan@cecm.sfu. ca or gjfee@cecm.sfu.ca. 

Below, we give some examples showing how EZ Face may be used. The left- 
aligned lines represent the input to EZ Face, while the centered lines represent the 
output of EZ Face. All computations are done with the precision of 50 digits. 

•^ Centre for Experimental an d Constructive Mathematics, Simon Fraser University, 
ittp ; //www . swox . com/gmp/ 
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Example 7.2. 

Pi-6/z(6) 

945 . 00000000000000000000000000000000000000000000000 

Example 7.3. 

lindep([z(4,l,3), z(5,3), z(8), z(5)*z(3), z(3) -2*z(2)] ) 

36. , 36. , -71. , 90. , -18. 

Example 7.4. 

lindep([ z(3), Pi-2*log(2), zp(2,2,l), zp(2,3) ]) 

12. , -1. , -12. , -12. 



Example 7^ is a simpl e in stance of Euler's formula for C(2n). Example 7^ 
the discovery of equation (|3.l[). Example 7.4 confirms formula ([7.4|). 



8. Evaluations for Unit Euler Sums 

As usual, the Holder conjugates p and q denote real numbers satisfying 1/p + 
1/q = 1, and p > 1 or p < —1 for convergence. Our first result is an easy conse- 
quence of the binomial theorem. 

Theorem 8.1. The generating function equality 

oo 
n=l 

holds. 



Proof. By definition ( p.3| ) of fj,, 



m— 1 



n=l m=l ^ j = l ^ -' 



D 



Corollary 1. 

M(M") = (logg)"/n!, 0<neZ. 

Remarks 8.2. Of course, when n = 0, we need to invoke the usual empty product 
convention to properly interpret /i({}) = 1- Since the mapping p i—^ 1 — p induces 
the mapping q i-^ 1/q under the Holder correspondence, duality (^) takes the par- 
ticularly appealing form /i({p}") = (— l)"/i({l— p}") in this context. In particular, 
p — —1 and (5-duahty (|6.8D, ( |6.10 ) gives 

-^({1}") - (-1)X{-1}") = (log2)7n!, < n e Z, 
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i.e. 

V TT - = V TT ^ ' 

fl,...,!/^ — Ij — 1 -^ ' l^i,...,Z^„ — Ij — 1 -^ 



(1°^2) 0<.ez, 



71! 

which can be viewed as an iterated sum extension of the weU-known resuh 



oo , oo / ^\v+\ 



typically obtained by comparing the Maclaurin series for log(l + x) when x = ~\ 
and a; = 1. 

We now prove a few results for unit Euler sums that were left as open conjectures 
in W. It will be convenient to employ the following notation: 

(8.1) A.:=Li.(i)=<5(r)=f;^, /^, :^ ii^, Z, := {-\)\{r). 

k=\ 

Theorem 8.3. For all positive integers m, 

771 

/i({-l}"\ 1) - (-1)™+! Y, Ak+lPm-k - Zm+l. 
k=0 



Proof. From the case (7^) of the Holder convolution, we have 



777+1 



<5(2, {l}'"-i) = C(™ + 1) - E 5{r)5{{iy^+'-n. 



Now multiply both sides by (— 1)™ and apply the case (3.11) of (5-duality. D 



Remarks 8.4. Theorem 8.3 appeared as the conjectured formula (67) in H, and 



is valid for all nonnegative integers m if the divergent m = Q case is interpreted 
appropriately. The equivalent generating function identity is 



' — ' Pi 

1 Jo 



t 



-dt 



_lo-2 V f ^— - - V V ^"'"''"^ 



.x + n nl ^-^ x + n 

n—l ^ ' 77—1 



correcting the misprinted sign in formula (21) of 



The asymmetry which marrs Theorem p. 3| is recovered in the generalization (2.4), 
restated and proved below. 

Theorem 8.5. For all positive integers m and all nonnegative integers n, we have 

M({-irM,{-ir) = (-i)™+iE(''^^)^fc+"+i^"-^ 
(8.2) + {-ir+^Y.h^^\zk+r.^+iPu-k. 



where A^, Pr and Zr are as in (8.1). 
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Proof. Let to be a positive integer, and let n be a nonnegative integer. We have 

M({-ir,i,{-ir) = (-ir+"+i / c^'_>i Tc^!!! 

Jo Jq 

= (-1)™+"+! f uj1\uJi I iO^ 

/■I /•(l-y)/2 

_ I 1 ^m+n+l , ,m n 

Jo Jl/2 

r-1 
\m+n+l I m 



(-1)"+"+' / ^L>i(log(l + y))'7"!. 

"'0 



By duality, 

m!n!/i({-ir,l, {-!}") = m\ f {- logi2 - y)r uocu^ 

Jo 
/■I ry/2 

= m! / (-log(2-2;))"c.o / < 
Jo Jo 

= /' (- log(2 - y)r (log(l - 2//2))" dj//y. 
Jo 

Letting t = I — y/2 and forming the generating function, it follows that 

oo oo 
771—1 n— 

= T.T.^\ (-iog(2t)r(iogir-^ 

f-^ (2i)-!' (r - 1) 



1/2 1 ^ ^ 



■dt. 



Expanding 1/(1 — t) in powers of t and integrating term by term yields 



m— 1 n— 



(8-3) =2-^E(^^ r^l-Ei^^^+E 



^ k + X ~ V k — y J ^-^ k + X — y ^-^ k — y 

k=l ^ ^ "^ fc=l ^ fc=l " 



Since m > 1, we may ignore the terms in ( |8.3| ) which are independent of x. Thus 
formally, but with the divergences coming only from the terms independent of x 
and hence harmless, 

oo f. oo 

_2-y-^ — + 2-'^y ^ 

Z-^ h 4- T — II ^ — ' 



k + x — y ^-^ k + X — V 

fe=i ^ fc=i *^ 

oo oo oo oo 

r=0 h=l r=0 h=l 



where we have used the abbreviations in (8.1). It is now a routine matter to extract 
the coefficient of a;™y" to complete the proof. D 
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Remark 8.6. Theorem 8.5 is an extension of conjectured formula (68) of W, and is 
valid for all nonnegative integers m and n if the divergent m = case is interpreted 
appropriately. 

9. Other Integral Transformations 

In Section o, we proved the duality principle for multiple polylogarithms by using 
the integral transformation y f—^ 1—x. Similarly, in this section we prove additional 
results for multiple polylogarithms by using suitable transformations of variables 
in their integral representations. 

Theorem 9.1. Let n be a positive integer. Let 6i,... ,bk be arbitrary complex 
numbers, and let si, . . . , Sk be positive integers. Then 

\Vl,b^,...,bl) ^ \exbu... ,ekbk 

where the sum is over all n^ cyclotomic sequences 

ei,... ,efc G {l,e2"/",e4'^*/",... , e^'^*^"-^)/"} , 
and, as usual, s := si + S2 + ■ ■ ■ + Sk. 
Proof. Write the left-hand side as an iterated integral as in (|4 

\0l,02, ■ ■ ■ ,0,./ Jo j^-^ 

Now let y = x" at each level of integration. This sends ujq to hloq and, by partial 
fractions, 

n-l 
r=0 

The change of variable gives 

^1 k n— 1 






2TTir/n 
j=l r=0 



Now carefully expand the noncommutative product and reinterpret each resulting 
iterated integral as a A-function to complete the proof. D 

Example 9.2. When n = 2 and fc = 1, Theorem |9.l| asserts that 

l + (-l)" 



C(.s)-2-i^ 



n=l 



Thus, Theorem 9.1 can be viewed as a cyclotomic extension of the well-known "sum 
over signs" formula for the alternating zeta function: 

E -^ = (1 - 2'"^)C(^)' ^(^) > 0- 

n—1 

Next we prove two broad generalizations of formulae (24), (26) and (28) of 0. 
By a pair of Cat operators we mean nested concatenation (similarly as two ^ signs 
mean nested summation). 
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Theorem 9.3. Let Si,S2, ... , Sfc be nonnegative integers. Then 

^. 1 + Sfc, 1 + Sk-l, ..., l + si\^\p 
-L, i, ■■■1 ^ / 



fi ( Cat{-1} Cat{e,,,}] fj f[ e,^ 



where the sum is over all 2*^+*'^'' '"^'^ sequences of signs (cij), ™t/i each ei, 
{1, —1} /or all 1 < i < Sj, 1 < j < k, and Cat denotes string concatenation. 

Proof. Let 

1 + Sfc, 1 + Sfc-i, ... , 1 + si ^ '■ 



L:=A ^+^^'^^' ^+j^|-i' ;;;' ^^^^ = (-i)^ n-o^--i 



,=fe 

Now let us use duality, and then we let y = 2t/(l + 1) at each level of integration. 
We get 

nl fe 

UJl) 

Now let us carefully expand the noncommutative product. We get 



Jo .1 






where the sum is over all sign choices Eij G {1,-1}, 1 < i < Sj, 1 < j < fc, and 
where by #eij = a we mean the cardinality of the set {(i,j) \ Si_j = a}. 

Let us now interpret the iterated integrals as A-functions. In this case, they are 
all unit Euler ^-sums, as we defined in (2.3). Thus, 

k 



L=i-l)''J2i-l)*^-^-\-l)''+^J Cat{-l}CatK,}], 



where, as usual, s := si + S2 + • • • + Sfc. Now if r of the Eij equal +1, then s — r of 
them equal —1. Hence, 

L = E(-1)^'"^="V( Cat{-1} c|tK,}) . 

Finally, (— l)#^i.j=^i is the same as the product over all the signs e^.j, and this 
latter observation completes the proof of Theorem 9.3. D 



Theorem 9.2 generalizes several identities conjectured in 0]. For example, we 
get the c onje cture (28) of [Q if we put s„+i = m, s„ — s„_i = ... = si = in 
Theorem 9^. Furthermore, (24) of Q is the case Sm+n+i — Sm+n — ■■■ — Sn+2 — 0, 
s„+i = I, Sn — Sn-i — ... — Si — 0, and (26) of is a special case of Theorem 3^ 
as well. Thus every multiple polylogarithm with all alternations (or, equivalently, 
every Euler sum with first position alternating and all the others non-alternating) 
is a signed sum over unit Euler sums. The representation of the sign coefficients 
used in Theorem 9.3 is much simpler than the cumbersome form of (28) in 0. 

Below we present a dual to Theorem |9.3| , which gives any unit Euler /i- value 
in terms of A- values with all alternations (equivalently, Euler sums with only first 
position alternating): 
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Theorem 9.4. Let si,S2, ■■■, Sk be nonnegative integers. Then 



J Cat{-l}{ip- J =Y.H Cat Cat {i,,,-} 



where the sum is over all 2''^^^'^^^ ^"'^ positive integer compositions 

tl,j + ^2j + ••• + tqj,o =5^+1) 1 < <l3 < S] +1, 1 < j < k. 

Proof. Let 

M:=M(^Cat{-l}{ir'=-^ =(-l)'=5(^Cat{l + s,}^ =^ H^o^c^s- 

Again, let us make the change of variable y = 2t/(l + t) at each level. Then 

.1 fe 



"'0 ,_-, 



'0 ,^1 

Again, let us carefully expand the nonconimutative product. Wc get 









i-^-l), 



where this time, the sum is over all Sij G {0, —1} with I < i < Sj, 1 < j < k. Note 
that each lu^i in the integrand contributes —1 to the sign and +1 to the depth. 
Since 

pi 
(— 1) P / weight-length string = A(depth-length string), 

it follows that M is a sum of A- values with all +1 coefhcients. That is, 

where the sum is over all vectors 

tj = (tij, . . . , tg^j), l<qj <l + Sj, 

and such that 

^tl,j = 1 + Sj, 1 < j < fc. 

i=l 

In other words, the sum is over all 2* independent positive integer compositions (in 
the technical sense of combinatorics) of the numbers I + Sj, I < j < k. D 

10. Functional Equations 

One fruitful strategy for proving identities involving special values of polyloga- 
rithms is to prove more general (functional, differential) identities and instantiate 
them at appropriate argument values. In the last two sections of this paper we 
present examples of such proofs. 
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Lemma 10.1. Let < a; < 1 and let 

Then 

(10.1) ,7(-.) ^ -J(.) + i,/(x') + ,/( Ji^) - h'{j^) ■ 



Proof. If £(x) and _R(x) denote the left-hand and the right-hand sides of ( 10.1 ), 
respectively, then by elementary manipulations (under the assumption < x < 1) 
we can show that dL/dx = dR/dx. The easy observation L{0) = i?(0) = then 
completes the proof. D 



Remarks 10.2. The identity (10.1) can be discovered and proved using a computer. 
Once the "ingredients" (the J-terms) of the identity are chosen, the constant coef- 
ficients at them can be determined by evaluating the J-terms at a sufficiently arbi- 
trary value of X e]0, 1[ and using an integer relation algorithm |10[ . Once the iden- 
tity is discovered, the main part of the proof (namely showing that dL/dx — dR/dx) 
can be accomplished in a computer algebra system (e.g., using the simplify () 
command of Maple) . 

Theorem 10.3. We have 

(10.2) A(2-,l-) = C(2,l)/8. 

Proof. Using notation of Lemma |l0.l| let us observe that 

Tii>n2>0 ^ 

Plugging in a; = 1 and applying (10.l|) now completes the proof. D 



Remarks 10.4. Theorem 10.3 is the n = 1 case of the conjectured identity (23) of 
|7J, namely 

(10.3) A(2-,l-,2,l,...)=8-"C({2,ir), 

" v ' 

2n 

for which we have overwhelming numerical evidence. This evidence also suggests 
that ( 10.3 ) with n > 1 seems to be the only case when two Euler sums that do 
not evaluate (in the sense of the definition in Section g|) have a rational quotient, 
different from 1. (See also Section p^.) 



11. Differential Equations and Hypergeometric Series 
Here, it is better to work with 

,Sk), 



,Si]X) 



i(sl,. 


■■ ,Sk]x) 


:= ^i/xisi, ■ ■ • 


since then we have 






l"-^''- 


■ ,si]x) = 


---Li-l + Sk, 


if Sfc > 2; while for s^ = 1, 






d 

Tx"-^''--- 


,si;x) = 


T L{sk^i, 
1 — x 



,si;x). 
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With the initial conditions 

L{sk, . . . , Si; 0) = 0, A: > 1, and L{{}; x) :— 1, 
the differential equations above determine the L-functions uniquely. 

11.1. Periodic Polylogarithms. If s := {si,S2,--- ,Sk) and s := J2j^j^ then 
every periodic polylogarithm L[{sY) has an ordinary generating function 



Ls{x,t):=Y,L{{sY-x)e 



r=0 

which satisfies an algebraic ordinary differential equation in x. In the simplest 
case, fc = 1, s reduces to the scalar s, and the differential equation for the ordinary 
generating function is Ds ~ f — 0, where 

The series solution is a generalized hypergeometric function 

r-l 



Ux,t) = i + E-^JiII 1 + 71 

r=l 7 = 1 ^ -^ ^ 



= s-Ts-l 



1,1,... ,1 

where w — e^'''^ , a primitive sth root of —1. 

11.2. Proof of Zagier's Conjecture. Let 2-^1(0, 6; c; a;) denote the Gaussian hy- 
pergeometric function. Then: 

Theorem 11.1. 



(11.1) EL({3,ir;2:)<4" 



n=0 

= 2F1 (ii(l + i), -it(l + i); 1; x) 2F1 (it(l - i), -it(l - i)- 1; x) . 
Proof. Both sides of the putative identity start 

^ + ¥^ +18^ +^53^^ +••• 
and are annihilated by the differential operator 

Once discovered, this can be checked in Mathematica or Maple. D 

Corollary 2. (Zagier's Conjecture) [|69| For all nonnegative integers n, 
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Proof. Gauss's 2^1 summation theorem gives 

1 sin(7ra) 

2Fi(a, -a; 1; 1) = -— — — — — - = . 

1 (1 — a)i (1 + a) na 

Hence, setting x = 1 in the generating function ( |l 1 . l| ) , we have 



00 

= 2F1 (it(l + i), -it(l + i); 1; 1) 2F1 (it(l - i), -it(l - i); 1; l) 

2sin(i(l + i)7rt)sin(i(l-i)7rt) 

" ^^ 

cosh(7ri) — cos(7ri) 

" ^ (4n + 2)!' 



n 



Remark 11.2. The proof is Zagier's modification of Broadhurst's, based on the ex- 
tensive empirical work begun in [[7|. 

11.3. Generalizations of Zagier's Conjecture. In H we give an alternative 
(combinatorial) proof of Zagier's conjecture, based on combinatorial manipulations 



of the iterated integral representations of MZVs (see Sections 4.2 and 5.4). Using 
the same technique, we prove in Ig] the "Zagier dressed with 2" identity: 

_-4n+2 

where s runs over all 2n + l possible insertions of the number 2 in the string {3, 1}". 
Still, ( 11. 2| ) is just the beginning of a large family of conjectured identities that we 



discuss in J8| 

12. Open Conjectures 

The reader has probably noticed that many formulae proved in this paper were 
conjectured in H. For the sake of completeness, we now list formulae from |R] that 
are stih open: (18), (23), (25), (27), (29), (44), and (70)-(74). It is possible that 
some of these conjectures can be proved using techniques of the present paper. 
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